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Recent experiments with superconducting qubits are motivated by the goal of 
fabricating a quantum computer, but at the same time they illuminate the 
more fundamental aspects of quantum mechanics. In this paper we analyze 
the physics of switching current measurements from the point of view of 
macroscopic quantum mechanics. 
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1. INTRODUCTION 

Since its very inception, quantum mechanics has defied our classical 
intuition. Quantum-mechanical correlations established between parts of a 
system during interaction are of a different nature than the classical onesP^ 
Much work has been put recently into harnessing the power of these corre- 
lations for performing computational tasks which are very difficult to imple- 
ment on classical computers. Superconducting qubits based on the Joseph- 
son effect have been proposecPas the elements of future quantum computers, 
based on the previously demonstrated macroscopic quantum coherence ef- 
fects in charge and flux devices P 

A number of superconducting qubits are currently under close experi- 
mental investigation, such as charge qubits^l, phase qubits^, flux qubits^, and 
a mixed charge-flux version called Quantronium. 7 Quantronium has a very 
large decoherence time (more than 500ns) , and it will be often referred to in 
the last part of this paper. To measure the qubit, a now standard technique 8 
is to monitor the switching probability of a large read-out junction or dc- 
SQUID to which the qubit is coupled. The switching probability depends 
on the state of the qubit, therefore a change in this probability at the same 
constant bias current indicates a different qubit state. 

If and when a quantum computer can be operated is an open question 
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that depends on the progress in reducing decoherence; even if if will be con- 
structed, a quantum computer will be able to address only a limited number 
of niche problems - such as factorization and database search - better than its 
classical counterpart. At the same time, the physics and technology behind 
present-day qubits is sound, and one may wonder whether these systems will 
have other potentially interesting applications besides quantum computing. 
Already technologies based on the properties of entangled light are at a ma- 
ture stage, leading to industrial applications, e.g. quantum cryptography. 
At the same time, fundamental research such as testing quantum mechan- 
ics at the macroscopic level is an important topic envisioned decades ago^, 
with progress in this direction enjoying now a firm experimental basis. The 
most spectacular test would be a clear experimental proof of violation of 
Bell's inequalities, which from the experimental point of view looks like a 
formidable task ahead of us requiring longer two-qubit decoherence times 
and read-out systems with a higher visibility than what is currently avail- 
able. Also, in a more general sense, quantum computing can be regarded as 
a test of quantum mechanics. 



2. THE JOSEPHSON EFFECT 

An underdamped Josephson junction^ can be described by the Hamil- 
tonian 

Q 2 

H = — - E j cos j - IcpoJ, (1) 

where I is the value of the bias current, 0o = ^o/2tt = h/2e = 3.295 x 1CP 8 
Gem 2 is the "barred" flux quanta, 7 is the phase difference of the supercon- 
ducting order parameter across the junction, Ej = (polco is the Josephson 
energy (J c q is the critical current), and C is the capacitance of the junction. 
This last electrostatic parameter is conveniently characterized by the energy 
associated with charging the junction capacitor with a single Cooper pair, 
Ec = (2e) 2 /2C. The dynamics of the junction is formally equivalent to sys- 
tems such as a particle (of "mass" C) in a washboard potential (see Fig. [TJ, 
or a gravitational pendulum under a constant torque. The degrees of free- 
dom of the electromagnetic environment in which the junction is embedded 
play an important role, but they will not be discussed in this paper. 

The charge Q and the magnetic flux <j)Qj (or equivalently the number of 
Cooper pairs n = Q/2e and the relative phase 7) are canonically conjugate 
variables. Hamilton's equations of motion give 

#07 _ Q , s 

dt ~ c { ' 
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Fig. 1. Energy levels associated with the washboard potential describing the 
dynamics of a curent-biased Josepshon junction (schematic presented to the 
right). 



^ = I-Ej^ S m r (3) 

The first equation Eq. is in fact Josephson's famous formula: a 
finite voltage across a junction is related to a change in phase. Equation 
([3]) is Kirchhoff's law for currents. The connection between this equation 
and Farady's law can be understood if one considers a superconducting ring 
interrupted by a Josephson junction. If we can neglect the inductance of the 
ring itself, then the quantity (ftoj equals the magnetic flux through the ring: 
according to Faraday's law, a change in this flux produces indeed a voltage 
(ftod'j/dt which equals the voltage across the junction Q/C. The second 
equation Eq. ([3]) expresses one of Kirchhoff's laws (conservation of charge in 
a node of an electrical circuit). The current Ij = Ejcf>Q sin 7 = / c osin7 is 
the current due to the Josephson effect. This current can be used to define 
a Josephson nonlinear inductance V = Ljdlj/dt, and using Eq. (|2|) we find 
that 

hi j cos 7 

Thus, a Josephson junction can be regarded as an LC-oscillator consisting of 
a flux-dependent nonlinear inductor in parallel with a capacitor. The inverse 
of the Josephson inductance Lj 1 measures the curvature of the washboard 
potential energy U{pf) = Ej(l — cos 7) + I(j)oj at any point 7 - in the same 
way as the inverse of capacitance is the second derivative of the charging 
energy Q 2 /2C. 
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The quantization of this system follows the usual recipe in quantum me- 
chanics: quantum-mechanical effects become important if the temperature 
is low enough (lower than the energy level separation) , which turns out to be 
the case for say typical Aluminum or Niobium junctions thermally anchored 
to the mixing chamber of a dilution refrigerator and with carefully filtered 
biased lines. The preferred variables to describe a quantized junction are 
the non-commuting charge and flux operators, [<j>oi, Q] = (or [7, n] = i). 
It can be readily checked that the Heisenberg equations of motion are for- 
mally identical (all variables being now understood as operators) with the 
Josephson-Kirchhoff relations Eq. ([2]) and Eq. ([3]). 

The physics of Josephson pendulums is not unique to metallic junctions: 
the pendulum Hamiltonian can be used in any situation in which two su- 
perfluids are connected by a weak junction. Consider for instance the case 
of bosonic atoms confined in a two-well potential, a situation which can be 
realized with alkali atoms trapped in optical lattices.^ The Hamiltonian of 
this system is a two-site Bose-Hubbard model with two components: an 
intrawell interaction energy and a tunneling part, 

H = -(a + a + aa + b + b + bb) - t(a + b + b+a), (5) 

where the annihilation operators a, b refer to the two wells. In this problem, 
the total number of particles N = a + a + b + b is a constant of motion. Clearly 
what happens in this system is that the dynamics can be described as an 
oscillation of the relative number of particles between the wells. In the 
regime w <C tN <C N 2 w (sometimes^ called the Josephson regime of a 
two-well Bose-Einstein condensate) one can introduce the variables 

a + a - b + b ^ 2 , ~ 
n = , e* 7 = —a+b, (6) 

we find that, up to constant terms, the Hamiltonian Eq. ([5]) assumes the 
form of a pendulum (without the bias term), 

H = -tN cos j + wn 2 . (7) 

In the presence of a bias current (or, in the case of alkali atoms, if the 
optical lattice is tilted in the gravitational field) the washboard potential 
[7(7) = — Ej cos 7— /0o7 allows a particle to tunnel in the semi-continuum of 
energy states immediately outside the well. How do we understand quantum- 
mechanically this process? Surely, this is an ubiquitous phenomenon which 
occurs for instance from nuclear a-decay to excited atoms emitting photons 
in vacuum: the physics is governed by the classical exponential law, which 
states that the fraction of undecayed particles after a time t is exp(— Ft), 
where T is the decay rate. 
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3. IRREVERSIBILITY 

For a Josephson junction, the decay is essentially a tunneling process 
between the state localized in one of the wells on one hand (which we call 
l^o(7)))) an d the states outside the well (Fig. [[]); the generic form for the 
macroscopic wavefunction at any time would then be^ 



|tf(t)) = e - r */2 e — ot|v£, o) + ^ _ e-^\^ out {t)). (8) 

Here the states ^out are approximately orthogonal to each other and 
normalized to 1. The decay amplitude exp(— rt/2), which, when squared, 
gives the correct classical decay probability law, is not straightforward to 
understand. Quantum-mechanically, the factor in front of a mode should 
be an exponential of an imaginary number (a product of energy and time 
divided by h), much like e~ tuJot . Also, the Schrodinger equation is time- 
reversible, while Eq. ([8]) is clearly irreversible. 

To understand where irreversibility results from, let us consider a simple 
modeP^in which the state |^o) can tunnel into the continuum of states {ipe} 
outside the barrier; the Hamiltonian for this model is 



H = huJo\^o)(^o\+ J defa# e )(^ e |+ J de [k(co , e)|* )(^e| + k*(u , e)|^ e )(* | 

(9) 

and the wavefunction can be expanded as 

= a(t)e-^ ot \* ) + [ deb(e, t)e" ie ' |t/> £ ). (10) 



This problem can be solved-^! using the Laplace transform and the clas- 
sical decay law results immediately. Here we would like to gain a better 
intuitive understanding of the mechanism that leads to the appearance of a 
decay from the otherwise reversible Schrodinger equation. Let us start by 
considering the Hamiltonian Eq. ([9]), this time written for a discrete set of 
out-of-the-well states {ip m }, with k(uo,e) — > &(o, m )) 

H = ?iUJo\^o}{^o\+J2 he m\^m)(lpm\+J2 [ fe (0,m) l*0> (^m| + fc( 0>TO ) l^m) (^0 1 • 
mm 

(11) 

The corresponding wavefunction expansion 

|*(t)> = a(i)e-^*|* ) + E M*)e~ iem Vm>, (12) 

m 

results, when used with the Schrodinger equation ihd\^{t))/dt = H\^/(t)), 
in an integrodifferential equation for a(t), 

da(t) 1 * 



dt h 



\ J ^T^ 0- ^"' ^)! 2 "^)- (13) 
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We now assume that the density of states (let us call it jV) and the 
tunneling amplitude fc(o,m) are slowly varying around the frequency Uq, near 
which the sum in the last part of the expression Eq. (113[> above is non-zero. 
This yields 



ps 2TThM(io )\k {m \ 2 5(t - t'). 
Inserting this result into Eq. (|13p we obtain 



^1 = -la(t), (14) 
dt 2 w ' v y 

where T = (27r /h)Af(uJo) |fy ,o) I 2 * s ^ ne decay rate. This set of approximations 
(sometimes called Weisskopf-Wigner theory-D) leads to the same formula for 
the decay rate as given by the Fermi golden rule (first order time-dependent 
perturbation theory) for transitions between a discrete level and a continuum 
of states. 

Let us now give a numerical illustration for this procedure. We take 
fyo,m) = fyo,o) energy-independent and we solve numerically the Schrodinger 
equation (with ft = 1 and &(o,o) = lj as a result time has no units). We 
consider 50 states such that e m — u>o = me(—l) m (in other words the states 
{IV'm)} are equally spaced around the resonance energy level set by ujq). Here 
e is the energy separation between the states (the inverse of the density of 
states). Clearly, keeping the number of states constant and going to lower 
and lower values of the energy separation e, the states tend to merge into 
a single level, resonant with |^o)- I n this 

CctSG, ELS expected, the particle 
oscillates between the intrawell state |^o) an d the level outside the well; in 
our simulation, we see this as oscillations between and 1 of the probability 
|a(i)| 2 . For larger values of e, the states began to separate and the dynamics 
is different. In Fig. [2] we present the probability |a(t) 2 | for e = 0.2. The 
inset shows the approximately exponential decay of the probability to find 
the system in the initial state |^o)- At larger time scales, some peaks ap- 
pear (corresponding to the particle returning to the well) due to accidental 
destructive interference of the amplitudes b m (t) corresponding to the states 
outside the well. This is a coincidental effect resembling the revival and 
collapses of the wavefunction in quantum optics; it can be removed either 
by increasing the number of states of by making the energy separation be- 
tween them not so uniform. In Fig. [3] we show the results of the second 
strategy, again with 50 states, fc(o,o) = 1> an d e = 0-2, but this time with 
e m — ujq = 50[RAND]e(— l) m , where [RAND] is a random fraction between 
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and 1. We notice that the possibility of large accidental constructive in- 
terference in the well is reduced, and also that the shape of the decay is very 
similar to that of the inset of Fig. [2] (we have checked this statement for 
various other values of [RAND]). 



a(t)\ 2 




5 10 15 20 



Fig. 2. The probability of finding the particle in the well for a model with 
constant spacing between the energy levels outside the well at large times. 
Inset: short-times detail showing the decay. 

We have then demonstrated that the Schodinger equation, albeit re- 
versible, can lead to an irreversible, decay-type evolution when a large num- 
ber of states is involved. This effect is related to the very small chance of 
having a constructive interference that would reconstruct the initial wave- 
function. 

4. THE QUANTUM MEASUREMENT PROBLEM 

Let us now consider the measurement problem for a typical supercon- 
ducting qubit coupled to a read-out system based on switching probabilities. 
The case of QuantroniumP is instructive in this sense. 

Here, a split Cooper pair box with Josephson energy Ej/2 per junction 
is coupled to a larger junction Ejq of electric capacitance Cq. To maximize 
the decoherence time, the qubit is kept at the optimal point for which the 
charging energy of the states \n = 0} and \n = 1) is degenerate (n counts 
the excess number of Cooper pairs on the island). The qubit states are of 
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Fig. 3. The probability of finding the particle in the well for a model with 
random spacing between the energy levels outside the well at large times. 
Inset: detail at short times displays almost the same decay law as in the 
previous figure. Note: the inset and the large graph have different values of 
the random parameter [RAND]. 



the Stern-Gerlach spin-1/2 type, given approximately by |+) 
and |— ) = 775 (|0) ~~ i n this basis the Hamiltonian is 

1 7 - Q 2 

H = --Eja z cos - - Ej cos 7 - J0 O 7 + 

The macroscopic wavefunction is then spin-dependent 

|*(7,t))=*+(7,t)|+) + *-(7,t)|-), 
where each component evolves according to 



yo)+ii» 



V2 



2q2 



h 2 d 



(15) 



(16) 



(17) 



with a spin-dependent wasboard potential U±( r y) = —Ejocosj + I(po7 T 
(Ej/2) cos(7/2). When the bias current is raised adiabatically close to the 
critical current Ejq^q 1 of the large junction, this leads to two tunneling 
rates T± which can be calculated from U±('y) by standard WKB methods. 
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Fig. 4. A Quantronium schematic. A Cooper pair box is connected by two 
Josepshon junctions to a larger current-biased read-out junction. 



We have showrP^ that for a measurement sequence in Quantronium, if 
we prepare the qubit in the state a|+)+/?|— ), then the macroscopic quantum 
state of the whole system evolves during the measurement sequence as 



o 



P 



+ + 



(-)/ 



(18) 



Here the normalization of ^f ut is to 1 - not the same as in the previous 
wor ]pl] w h ere ^ k ee p the equations shorter we have normalized ^ ou t to [1 — 
exp(— Tt)]. This equation leads immediately to a formula for the switching 
probability during an interval r, 



P ou \t) = 1 



|a| 2 e 



|0[ 5 



(19) 



Equation (I19h above describes precisely the oscillations seen in the ex- 
periment, with a and dependent on the duration of the microwave pulse 
(as sine and cosine functions, with Rabi frequency set by the microwave 
intensity^) . 

What type of measurement does Eq. (I18p describe? Let us consider 
the case of ideal visibility (not yet achieved experimentally), namely the 
situation in which the states |— } and |+) can be distinguished 100% by the 
switching or non-switching of the large junction: this means exp(— T_r) = 
0, and exp(— r + r) = 1. The switching probability therefore simplifies to 
P out ( T ) = \f3\ 2 = 1 - |a[ 2 e- r + T , and the state Eq. ([18]) becomes 



o c : 



^ +) (7)l+)+/?^(7,t)l 



(-)/ 



(20) 
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This is precisely the form of a von Neumann picture of quantum measure- 
ment. Initially the qubit is prepared in a superposition of |— ) and |+), and 
the read-out junction is in the ground state, and during the "interaction" 
between the qubit and the junction (due to lowering the potential barrier 
by increasing the bias current), the two subsystems become entangled. The 
next step is the collapse of the wavefunction onto one of the components 
^ / q + ^(7)|+) or ^o U l (7, t)\— ), with corresponding probabilities \a\ 2 and \/3\ 2 . 

One thus immediately sees that Josephson-based superconducting quan- 
tum circuits are a perfect tool to study the conceptual foundations of quan- 
tum mechanics. A first observation is that the collapse of the wavefunction 
is a necessary ingredient in the fabric of the quantum formalism, albeit a 
rather unnatural one from a conceptual point of view. There is no way to 
get around this postulate by involving more and more subsystems between 
the qubit and the experimentalist, as it should be clear from the case of 
Quantronium, where the large junction does precisely the "interpolation" 
between the qubit and the classical world. The same argument goes also if 
one involves the degrees of freedom of (typically) the electromagnetic envi- 
ronment: to obtain the "real", classical probabilities, we are simply bound 
to collapse a more complicated wavefunction. The conceptual problem here 
is that the choice of states onto which the collapse happens is external to 
the theory itself. In our case, we know that the experimentalist either sees 
or does not see a switching event (by the appearance of a relatively large 
spike in the voltage across the junction, easily measurable by a voltmeter), 
therefore we know that for example *$> ut will be correlated with the classical 
movement of the (macroscopic) voltmeter indicator (pointer). But this is 
not a consequence of the theory: it is rather an external ingredient. Indeed, 
quantum theory allows us equally well to do the collapse on a superposition 
of ^0 and ^out- This observation is at the core of the issues associated with 
the Schrodinger cat measurement paradoxes: the punch line of this famous 
argument is that, if we insist that the quantum- mechanical wavefunction is 
not just a mathematical instrument to describe the outcomes of a given ex- 
periment but is "real" in the same sense as the waves on a lake are, then we 
must admit that superpositions of macroscopic objects could exist. The ex- 
perimental test of this possibility is of great importance for the foundations 
of quantum theory, and Josephson circuits are now becoming an essential 
instrument in the toolbox of the future quantum mechanic. 
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5. CONCLUSIONS 

We have investigated the physics of switching in Josephson junctions 
and we have described how this type of irreversible behavior emerges from 
the Schrodinger dynamics on a large number of states. For qub it-junction 
systems, we have analyzed the measurement process and we have shown 
that it is of von Neumann type. For superconducting circuits, as with all 
other quantum-mechanical systems, standard quantum theory allows us to 
describe quantitatively how Schodinger cats die, but offer no insight into the 
problem of how they choose the states onto which to collapse. 
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